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Abstract 

In the context of Calabi-Yau string models we explore the origin of char- 
acteristic pattern of quark-lepton masses and the CKM matrix. The discrete 
i?-symmetry Zk x Z2 is introduced and the Z2 is assigned to the i?-parity. 
The gauge symmetry at the string scale, SU(6) x SU(2)r, is broken into the 
standard model gauge group at a very large intermediate energy scale. At 
energies below the intermediate scale down-type quarks and also leptons are 
mixed with unobserved heavy states, respectively. On the other hand, there 
are no such mixings for up-type quarks. Due to the large mixings between light 
states and heavy ones we can derive phenomenologically viable fermion mass 
hierarchies and the CKM matrix. Mass spectra for intermediate-scale matter 
beyond the MSSM are also determined. Within this framework proton lifetime 
is long enough to be consistent with experimental data. As for the string scale 
unification of gauge couplings, however, consistent solutions are not yet found. 



1 Introduction 



In order to make sure of the reality of string theory it is important to explore 
how string theory determines low-energy parameters which are free parameters in 
the standard model and in the minimal supersymmetric standard model (MSSM). 
Among many issues of low-energy parameters, the characteristic pattern of quark- 
lepton masses and mixing angles has long been a challenging problem to explain its 
origin. The observed masses of quarks and leptons have the hierarchical pattern 

(i) m(lst gen.) <C m(2nd gen.) <C m(3rd gen.) 
and also the ratios among quark masses are in line as 

(ii) m u /m d < m c /m s < m t /m b . 

Up to now several possibilities of explaining these features have been studied by many 
authors A possibility is that all the observed pattern of fermion masses are 

attributable to the boundary condition, i.e. to the hierarchical structure of Yukawa 
couplings themselves at a very large scale. However, when we take GUT-type mod- 
els, it is difficult to find a satisfactory solution in which property (ii) comes into line 
with a simple unification of Yukawa couplings. In this paper we explore a somewhat 
distinct possibility. In the context of Calabi-Yau string models with Kac-Moody 
level-one we propose a new type of model which potentially generates the charac- 
teristic pattern of fermion masses and Cabibbo-Kobayashi-Maskawa(CKM) matrix. 
In the model property (i) is attributed to the texture of renormalizable and nonre- 
normalizable interactions restricted by some discrete symmetries at the string scale. 
This mechanism is similar to those proposed in Refs. PJ [0] H H] ■ On the other hand, 
property (ii) comes from large mixings among states observed at low energies and 
unobserved heavy ones. The mixings occur only for down-type quarks and for leptons 
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below the energy scale at which the gauge group is broken into the standard model 
gauge group G st = SU(3) C x SU(2) L x U(1) Y . 

The four-dimensional effective theory from Calabi-Yau string compactification is 
far more constrained than ordinary field theory. In the effective theory there are many 
peculiar features beyond the MSSM. First point is that the gauge group G, which is 
given via the flux breaking at the string scale Ms, is a subgroup of Eq and would be 
larger than the standard model gauge group G st . We will choose G = SU(6) x SU (2)r, 
under which doublet Higgs and color-triplet Higgs fields transform differently H. As 
we will see later, the gauge group G is spontaneously broken to G st in two steps at 
very large intermediate energy scales. Second point is that the massless sector of 
the Calabi-Yau string model contains extra particles beyond the MSSM. In string 
inspired models we typically have a number of generations and anti-generations. For 
illustration, if the gauge group G is E 6 , the massless chiral superfields apart from 
^-singlets consist of 

N f 27 + 5(27 + 27*), (1.1) 

where Nf means the family number at low energies. It should be noted that S sets 
of vector-like multiplets are included in the massless sector. In Calabi-Yau string 
compactification the generation structure of matter fields is closely linked to the 
topological structure of the compactified manifold. We will assume 6 = 1 for the 
sake of simplicity. Adjoint Higgs representations which are introduced in the tradi- 
tional GUT are not allowed at Kac-Moody level-one. In addition, particles beyond 
the MSSM are contained also in 27-representation of E$. Namely in 27 we have 
quark superfields Q = (U,D), U c , D c , lepton superfields L = (N,E), N c , E c , Higgs 
doublets H u , H d , color-triplet Higgses g, g c and an S'O (lO)-singlet S. When the 
gauge group G is broken into G s t, superfields D c and g c as well as L and Hd are 
indistinguishable from each other under G st . Hence there possibly appear mixings 
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between D c and g c and between L and H^. On the other hand, for up- type quarks 
there appear no such mixings. While up-type, down-type quarks and leptons share 
their interactions in common at the string scale Ms, D c -g c mixing and L-Hj mixing 
potentially turn mass pattern of down-type quarks and leptons out of that of up-type 
quarks at low energies. Further the mixings may be responsible for the CKM matrix. 
Third point of peculiar features beyond the MSSM is that superstring theory natu- 
rally provides the discrete symmetries which stem from symmetric structure of the 
compactified space. As shown in Gepner model |J, the discrete symmetry can be the 
.R-symmetry under which the components of a given superfield transform differently. 
Also, the discrete symmetry could be used as a horizontal symmetry. The discrete 
.R-symmetry strongly limits the renormalizable and nonrenormalizable interactions 
and then possibly controls parameters in the low-energy effective theory. Recently it 
has been argued that the discrete .R-symmetry controls energy scales of the symme- 
try breaking |TD[, the magnitude of Majorana masses of the right-handed neutrino 



j I !| and the stability of the weak-scale hierarchy [ E| . We will introduce the discrete 
.R-symmetry Zk x Z 2 at the string scale. 

In this paper main emphasis is placed on how both of the mixing mechanism 
mentioned above and the discrete symmetry bring about phenomenologically viable 
fermion mass pattern and the CKM matrix. This paper is organized as follows. In 
section 2 we introduce the discrete .R-symmetry, which puts stringent constraints 
on interactions in the superpotential. The Z 2 symmetry is chosen so as to be in 
accord with the so-called .R-parity in the MSSM. .R-parity is conserved over the 
whole energy range from the string scale to the electroweak scale. After arguing that 
the discrete .R-symmetry controls energy scales of the gauge symmetry breaking, we 
study particle spectra of vector-like multiplets in section 3. Since doublet Higgses 
and color-triplet Higgses belong to the different representations of the gauge group 
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G = SU (6) x SU (2)r, distinct particle spectra of these fields are derived without some 
fine-tuning of parameters. In section 4 mass matrices for colored chiral multiplets 
are presented. There appear mixings between D c and g c . Choosing appropriate 
assignments of discrete charges, we get large mixings between them. Due to the 
maximal mixing the mass pattern of down-type quarks differs from that of up-type 
quarks. The model generates not only the hierarchical pattern of quark masses but 
also the texture of the CKM matrix. In section 5 we discuss mixings between L and 
Hd and study spectra of leptons. The CKM matrix for leptons turns out to be a unit 
matrix. In the present framework we have several i?-parity even colored superfields 
which potentially mediate proton decay. In section 6 it is shown that proton lifetime 
is about 10 33 ~ 35 year. In section 7 we find that the gauge coupling unification is not 
successfully achieved as a consequence of spectra of extra intermediate-scale matter. 
In the final section we conclude with a brief summary of our results. In Appendix A it 
is shown that under an appropriate condition on the soft SUSY breaking parameters 
the gauge symmetry is broken at tree level. In Appendix B we show that if neutrino 
Majorana masses are sufficiently large compared with the soft sup ersymmetry( SUSY) 
breaking scale m 3 / 2 = O(lTeV), the scalar potential is minimized along the direction 
where i?-parity is conserved. 

2 Discrete i?-symmetry 

In order to guarantee the stability of the weak-scale hierarchy without fine-tuning, 
it is favorable that doublet Higgses and color-triplet Higgses reside in different irre- 
ducible representations of the string scale gauge group G. As the largest gauge group 
implementing such a situation is G = SU(6) x SU(2) R ||, in this paper we choose 
SU(6) x SU(2)r as an example of G. Chiral superfields (<&) in 27 representation of 
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Eq are decomposed into 

$(15,1) : Q,L,g,g c ,S, 

$(6*, 2) : U c ,D c ,N c ,E c ,H u ,H d . 

Although L and H d (D c and g c ) have the same quantum numbers under G s t, they 
belong to different irreducible representations of SU (6) x SU (2) R . The superpotential 
W is described in terms of 27 chiral superfields (<£>) and 27* ones (<&) as 

W = $ 3 + $ 3 + + $ 3 ($$)« + . . . ? (2.1) 

where m and n are positive integers and all the terms are characterized by the cou- 
plings of 0(1) in units of M s = O(10 18 GeV). The cubic term $ 3 is of the forms 

($(15, l)) 3 = QQg + Qg c L + g c gS, (2.2) 
$(15,1)($(6*,2)) 2 = QH d D c + QH U U C + LH d E c + LH U N C 

+SH u H d + gN c D c + gE c U c + g c U c D c . (2.3) 

We assume that the massless matter fields are composed of chiral multiplets 
(i = l,---,Nf = 3) and a set (5 = 1) of vector-like multiplets $ and Here 
we introduce the discrete i?-symmetry Z K x Z 2 as a stringy selection rule. As we 
will see below, large K is favorable for explaining the mass pattern of quarks and 
leptons. The Z 2 symmetry is taken so as to be in accord with the i?-parity in the 
MSSM. Therefore, hereafter the Z 2 symmetry is referred to as i?-parity. Supposing 
that ordinary quarks and leptons are included in chiral multiplets (i = 1,2,3), 
i?-parity of all $j (i = 1,2,3) are set to be odd. Since light Higgs scalars are even 
under i?-parity, light Higgs doublets are bound to reside in <3> and/or <3>. For this 
reason we assign even i?-parity to $o an d In Appendix B we show that once the 
.R-parity is conserved at the string scale, the i?-parity remains unbroken down to 
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the electroweak scale under appropriate conditions. Hence, through the spontaneous 
breaking of gauge symmetry gauge superfields are possibly mixed with the vector-like 
multiplets $o an d $ but not with the chiral multiplets (i — 1, 2, 3). Furthermore, 
no mixing occurs between the vector-like multiplets and the chiral multiplets. 

We use the Z% symmetry as a horizontal symmetry and construct our model 
incorporating the mechanism of Froggatt and Nielsen 0. The symmetry controls 
not only a large hierarchy of the energy scales of the symmetry breaking but also the 
texture of effective Yukawa couplings. We denote the -Z^-charges of chiral multiplets 
<3>i(15, 1) and $j(6*, 2) by a, and bi (i = 0, 1, 2, 3), respectively. In Table I, we tabulate 
the notations for Z^-charges and the assignment of i?-parity for each superfield. Note 
that the anticommuting number 9 has also a Z^ x ^-charge (— 1 ; — ). 



Table I 



3 Gauge hierarchy and the /i-term 

The discrete symmetry introduced above puts stringent constraints on both renor- 
malizable and nonrenormalizable interactions in the superpotential. To begin with, 
Z^-charges of vector-like multiplets are chosen such that both the nonrenormalizable 
interactions 

($ (15,l)$(15*,l)) Sfc (3.1) 

and 

($ (6*,2)$(6,2)) S (3.2) 

are allowed, where K — sk + 1 and s and k are even and odd integers larger than 
unity, respectively. This implies that sk(a Q + a) + 2 = s(b + b) + 2 = in modulus 
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K = sk + 1. Thus we impose 



a + a = 2 



b + b = 2k 



mod 



(3.3) 



It follows that the interactions 



s 



t^v ~ E( $ o(15, 1)$(15*, l))( s "^(*o(6*, 2)$(6, 2)) r 



(3.4) 



are allowed in Ms units. Due to i?-parity conservation the interactions containing 
even number of $j (i = 1, 2, 3) are also allowed but all of the interactions containing 
odd number of (i = 1, 2, 3) are forbidden. 

Incorporating the soft SUSY breaking terms together with the F- and D-terms, 
we get the scalar potential. Although dynamics of SUSY breaking is not presently 
known, we may parametrize the SUSY breaking by introducing the universal soft 
terms. The scale of SUSY breaking m 3 / 2 is supposed to be O(lTeV). Through the 
minimization of the scalar potential we are able to detemine a ground state, which is 
characterized by VEVs of <&q, $ and $j (i = 1, 2, 3). Under appropriate conditions on 
soft SUSY breaking parameters the gauge symmetry is spontaneously broken at tree 
level (see Appendix A). Further, if masses of G st - neutral and i?-parity odd superfields 
are sufficiently larger than m 3 / 2 , the scalar potential is minimized at vanishing 
for i = 1, 2, 3 (see Appendix B). On the other hand, $ and $ acquire nonzero VEVs 
along a D-flat direction, namely 



($0(15,1)) 



($(15*, 1 



)) 



~ M s x, 



(3.5) 



($o(6*,2 



)) 



($(6,2)) 



M s x' 



,k 



(3.6) 



up to phase factors |10H [p^Tfl , where 




(3.7) 
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Although for a large K the parameter x by itself is not a very small number, the large 
hierarchy occurs by raising the number to large powers. Hence, x becomes an efficient 
parameter in describing the hierarchical structure of the effective theory. Note that 
we have the inequalities 

M s > |($ (15,1)>| > |($ (6*,2)>| > ^m 3/2 M s . (3.8) 

Hereafter the fields $ (15, 1) and $(15*, 1) which develop non-zero VEVs are referred 
to as G si -neutral fields S and S, respectively. At the scale (So) = (S) ~ Msx the 
gauge symmetry SU(6) x SU(2)r is spontaneously broken to SU(4:)ps x SU(2)l x 
SU(2) R , where SU(A) PS stands for the Pati-Salam SU(A) @. Under the SU(4) PS x 
SU(2)l x SU(2)r the chiral superfields $(15, 1) and $(6*, 2) are decomposed as 

(15,1) = (4, 2,1) + (6, 1,1) + (1,1,1), (3.9) 
(6*, 2) = (4*, 1,2) + (1,2, 2), (3.10) 

where each matter field is assigned as 

$(4,2,1) : Q,L, 

$(6,1,1) : g,g c , 

$(1,1,1) : S, 

$(4*, 1,2) : U C ,D C ,N C ,E C , 

$(1,2,2) : H u ,H d . 

The subsequent symmetry breaking takes place via the non-zero VEVs ($o(6*, 2)) 
= ($(6, 2)) ~ Msx k . At this stage of the symmetry breaking there seem to be two 
possibilities depending on whether the non-zero VEV ($ (6*, 2))(($(6, 2))) is at- 
tributed to ($ (4*,1,2))(($(4,1,2))) or ($ (1, 2, 2))(($(1, 2, 2))). As will be seen 
soon later, we have the term (S S) p SH u H d in the superpotential, where p is a pos- 
itive integer detemined by the discrete symmetry Z^- Under an appropriate charge 
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assignment of matter fields we have p ~ sk — 2k. According to the presence of 
this superpotential term the large VEV ($(1, 2, 2)) is inconsistent with the (almost) 
F-flat condition. Consequently, the subsequent symmetry breaking occurs through 
($ (4*,1,2)) = ($(4,1,2)) ~ M s x k . Then we denote the fields $ (4*,1,2) and 
$(4, 1, 2) with the non-zero VEVs as Nfi and iV c , respectively. Thus the gauge sym- 
metry is spontaneously broken in two steps at the scales (So) and (Nq) as 

SU(6)xSU(2) R SU(4) PS x SU(2) L x SU(2) R (3.11) 

SU(3) C x SU(2) L x U(1) Y . (3.12) 

Since So, S, N£ and N° acquire VEVs along a D-flat direction, SUSY is maintained 
down to O(lTeV). 

At the first step of the symmetry breaking chiral superfields Qo, Lq, Q, L and 
(So — S)/y/2 are absorbed by gauge superfields. Through the subsequent symmetry 
breaking chiral superfields Uq, Eq, U°, E° and (Nq — N c )/\/2 are absorbed. On the 
other hand, for components (S + S)/a/2 and (Nq + N c )/y/2 the mass matrix is of 
the form 

0(x 2sk - 2 ) 0(x^ s -^ k - x ) \ 



\ 0(x^ 2s ^ k - 1 ) 0(x 2 ^ k ) 
in Ms units. This yields mass eigenvalues 



(3.13) 



0(m 3/2 ), 0(M s x 2(s -^ k ), (3.14) 



which correspond to the eigenstates 



j=(S + S) + O(x k - 1 )j= 2 (Nc + N c ), 
i l atc i ~Kf c \ i n(„k-l\ l 



(3.15) 



±(NS + AT) +O(x k ~ 1 )^(S + S), 

respectively [PTTf] . The discrete symmetry Z K is broken together with SU(6) x SU(2) R 
by the VEV (So), while the VEV allows the Z 2 -symmetry (referred to .R-parity con- 
servation) to remain unbroken all the way down to TeV. 

10 



In order to stabilize the weak-scale hierarchy we put an additional requirement 
that the interaction 

(S SY k - e S H u0 H d0 (3.16) 

is allowed with e = 0, 1 in the superpotential. We will shortly show that the \i- 
problem is solved by this setting e = 0, 1. This condition is translated into 

a + 2b = 2e mod K. (3.17) 



From Eqs. (|3.3f ) and ( 3-171) the superpotential of Higgs doublet in vector-like multiplets 
has the form 

+ (S S) sk ~ e SoH u0 H do . (3.18) 

When Sq and S develop the non-zero VEVs, the superpotential induces the mass 
matrix of H u0 , H d0 , H u and H d 

H u H d0 

H d ( 0{x 2 ^ k+2e ~ l ) 0(x 2 ^ k ) \ ( 3 ' 19 ) 
H u0 [ 0(x 2 ^ k ) 0(x 2sk - 2e+1 ) ) 

in Ms units, which leads to the mass eigenvalues 

0(M 5 x 2(s - 2)fc+2e " 1 ), 0(M s x 2sk - 2e+1 ) = 0(m 3/2 x 3 ~ 2e ). (3.20) 

Consequently, we have the //-term with m 3 / 2 > fi = 0(m 3 / 2 x 3 ~ 2e ) = O(10 2 ~ 3 GeV) 
for e = 0, 1 || [Q. Here, note that we take x ~ 0.7 with sk = 50 in a typical example 
given later. Light Higgs states are given by 

H u0 + 0(x 2k+1 ~ 2e )H d , H d0 + 0(x 2k+1 - 2e )H u . (3.21) 

The components of H d and H u in light Higgses are small. Generally speaking, in the 
superpotential Wh there exist additional terms which are obtained by replacing each 
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factor (SoS) k by a factor (NqN c ). However, as far as the mass matrices are concerned, 
these terms yield the same order of magnitude as in each entry of the above matrix 
because of the relations (\(S )\/M s ) k = \ (Nq)\/M s and k(a + a) = b + b. Since we 
do not address here the issue of CP-violation, all VEVs are assumed to be real for 
simplicity. Therefore, hereafter the nonrenormalizable terms are expressed in terms 
only of the powers of (SqS). Note that the product H u oHdo has a nonzero Z^-charge. 



In contrast with the present model, in a solution of the /i-problem proposed in Ref . [JT4 
the .R-charge of the product of light Higgses has to be zero. 

The remaining components in <3> and $, i.e. go, Dq and g, g c } D° are down- 
type color-triplet fields. In the present model the spectra of color-triplet Higgses are 
quite different from those of doublet Higgses. Mass matrix for these fields is given in 
section 6. 

4 Quark masses and the CKM matrix 

Next we turn to mass matrices for chiral multiplets (i = 1,2,3). Due to 
.R-parity conservation (i = 1,2,3) are not mixed with vector-like multiplets $o 
and $. The superpotential of up-type quarks which contributes to the mass matrix 
of up-type quarks, is given by 

W v ~ (SoSr-Q^H^ = 1,2,3), (4.1) 

where the exponents m^- are integers in the range < rriij < K = sk+1. Although the 
Zk symmetry allows the terms multiplied by (SqS) k , (S S) 2K , • • •, the contributions 
of these terms are negligibly small compared with the above ones. Therefore, it is 
sufficient for us to take only the terms with m^- < K. Recall that light Higgs doublets 
are almost H u q and H^o- Under the Z^-symmetry the exponent is determined 
by the condition 

2mij + Oi + bj + b + 2 = mod K. (4.2) 
12 



Instead of aj and bi (i = 1,2,3), hereafter we introduce new notations a, (3, 7, and 
S defined by 



a 2 — ai = a, b 2 - h = (3, a 3 - a 2 = 7, b 3 -b 2 = 5. 



(4.3) 



These parameters are supposed to be even integers to derive desirable mass pattern 
of quarks and 

< a < 5 < 7 < 0, 2(13 + 5) < K. (4.4) 

The above condition (|4.2j ) is rewritten as 

/ a + (3 + 7 + 5 a + 7 + 5 a + 7 



2mj,- = 2m 33 + 



/? + 7 + ^ 7 + <5 
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mod K, 



(4.5) 



where 2777,33 = — a 3 — b3 — bo — 2. The mass matrix of up-type quarks is described by 
a 3 x 3 matrix M with elements 

M ij = 0{x 2m ^) (4.6) 

multiplied by v u = (H u0 ). This equation is an order of magnitude relationship, so 
that each element will be multiplied by an 0(1) number. From Eq. Q4.5p the matrix M 
is generally asymmetric. Here we take an ansatz that only top-quark has a trilinear 
coupling. This means that 

m 33 = 0. (4.7) 
When we adopt appropriate unitary matrices V u and U u , the matrix 

V7 a l MU u (4.8) 

becomes diagonal. Explicitly, V u and U u are of the forms 



( 1 - 0(x 2a ) 

0(x a ) 
\ 0(x a+ ~<) 

1 - 0(x 213 ) 
0{x") 
Oix^ 5 ) 



0(x a ) 
1 - 0(x 2a ) 
0(x~f) 

0{x p ) 
1 - 0(x 2S ) 
0(x 5 ) 



0(x a+ ^) 

o(x~<) 

1 - 0(x 2 ~<) 

0(xP +s ) 
0(x s ) 
1 - 0(x 25 ) 



(4.9) 
(4.10) 
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The eigenvalues of M are 



0(x a+ ^ +s ), 0(x^ +6 ), 0(1), 



(4.11) 



which correspond to u-, c- and t-quarks, respectively. 

Under SU (6) x SU (2) R gauge symmetry down-type quarks and leptons share the 
nonrenormalizable terms in common with up-type quarks. Namely we get 



For down-type quarks, however, the mixings between g c and D c should be taken into 
account at energies below the scale (Nq). This is because we have two down-type 
S , [/(2) i -singlet colored fields in each 27 of E 6 . Then, hereafter we denote -R-parity 
odd gi and gl (i = 1, 2, 3) as D[ and D'f {i = 1, 2, 3), respectively. The superpotential 
of down-type colored fields is of the form 



W ~ (SoSrHQ^Hw + LtN'Huo + L t E]H d ^\. 



(4.12) 



W D ~ (SoSy-SoDlDf + (S S) m ^(N^D' t + H^D], 



(4.13) 



where the exponents % are determined by 



2z^ + di + cij + a + 2 = 



mod K 



(4.14) 




(4.15) 



Z tJ =0(x 2z -), 



(4.16) 



a mass matrix of down-type colored fields is written as 



jyc 




(4.17) 
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in Ms units below the scale (Nq), where pd = (Hdo)/Ms = Vd/Ms- This Md is a 6 x 6 
matrix and can be diagonalized by a bi-unitary transformation as 

Vd'MdUd- (4.18) 

Md shows mixings between D' c and D c , explicitly. This type of mixings does not 
occur for up-type quarks. From Eqs. ( 4.17|) and (|4.18| ) the matrix 



V^M d MlV d = V,- 1 ( Ad + B ^ *g ) V d (4.19) 

is diagonal, where 

A d = x 2 ZZ\ B d = x 2k MM\ 6 d = p d x- k . (4.20) 

In view of the smallness of the parameter e^, we use the perturbative method in 
solving the eigenvalue problem. It follows that the eigen equation is approximately 
separated into two pieces. For heavy states the eigen equation becomes 



det^Ad + Bd-^) 0. (1.21) 



Solving this equation of a variable rj, we obtain masses squared for three heavy 
states. The other three states are light and their masses are given by solving the 
eigen equation 

det(x- 2k (A-/ + B d ^-^j=0. (4.22) 

This equation is derived in e 2 d order of the perturbative expansion. The light states 
correspond to observed down-type quarks. If the mixing between D' c and D c is 
sizable, mass pattern of down-type quarks is possibly changed from that of up-type 
quarks. Thus in our model, property (ii) pointed out in section 1 for observed fermion 
masses is attributable to this mixing mechanism. 
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The 6x6 unitary matrices V d and Ud are 

V ~ ( Wd -e d (A d + B d )-'B d Vd\ , v 

Vd ~ [e d B d (A d + B d )- 1 W d V d J' (423) 

d " U fc MtW d (Ai°V 1/2 (x fc M)- 1 V d (Aj ) ) 1 /2 J ' l4 ' Z4j 
respectively. Here and V d are 3x3 unitary matrices which are determined such 
that the matrices 

W d \A d + B d )W d = Af, VJ W + B^Y'Vd = A$> (4.25) 

become diagonal. As a consequence we can expect to have a nontrivial CKM matrix 

V CKM = V~ l V d . (4.26) 

Note that V u is determined such that V~ l B d V u is diagonal. If the relation 

| det(Ai + B d )\~ \detA d \ > | det £ d | (4.27) 

is satisfied, the mixing is small and we have 

(A, 1 + B, 1 )- 1 - B d . (4.28) 

This implies that mass pattern of down-type quarks is the same as that of up-type 
quarks and that Vd — V u . In this case V CKM becomes almost a unit matrix. 

To get a phenomenologically viable solution, large mixings between D' c and D c 
are preferable. Thus we impose the maximal mixing in which (Aj 1 ),^ and (B d r )ij 
are the same order. The maximal mixing is realized under the condition 

2z 33 = jfe-l- a -7 + /3 + £. (4.29) 

Note that k is an odd integer. Under the above condition on Z33 the eigenvalues of 
A d + B d become 

0(:r 2(fe+a+7+5) ), 0(:r 2(fe+/3+5 - a) ), 0[x 2k ). (4.30) 
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It follows that extra down-type heavy quarks have their masses 

M sX k+a+1+ \ M s x k+0+5 ~ a , M s x k . (4.31) 

Main components of these eigenstates are Di-(0(l)Dl+0(l)D%), D 2 -D'^ and D^-D'f, 
respectively. On the other hand, down-type light quarks have their masses 

v d x a+f3 ^ +s , v d x^ +s , v d x~ a+l3+5 , (4.32) 

which correspond to observed d-, s- and 6-quarks. These eigenstates are approxi- 
mately Di-(0(1)D; + 0(1)Z>§), D2-D1 and D 3 -D' 3 C , respectively. It should be noted 
that we have very large Df-D'f mixings. The unitary matrix V d which diagonalizes 

A^ 1 + Bj 1 , is expressed as 

/ 1 - 0(x 2a ) 0(x a ) 0(x a+ ~<) \ 
V d = 0(x a ) 1 - 0{x 2a ) 0{x^) . (4.33) 
V 0(x a+ ^) O(x^) l-0(x 2 ~<) ) 

Corresponding elements of the matrices V u and V d are in the same order of magnitudes 
but their coefficients of the leading term in off-diagonal elements are different with 
each other because of the maximal mixing. Consequently, the CKM matrix is given 
by 



yCKM = V -1 V 



/ 1 - 0(x 2a ) 0(x a ) 0(x a+ ^ 



0(x a ) 1 - 0(x 2a ) 0{x-y) I . (4.34) 
V 0{x a+ ~<) 0(x 7 ) l-O(ar^) 

It is worth noting that large Df-Df mixings play an essential role in generating a 
nontrivial CKM matrix. An early attempt of explaining the CKM matrix via D^-D'f 
mixings has been made in Ref.|]l5[, in which a SUSY 50(10) model was taken. 



Confronting the CKM matrix obtained here with the observed one, it is feasible 
for us to take 

(t = 1.0xtc, 7 = 2.flxt« (4.35) 

with x w = A = sin 9c, where 9c is the Cabbibo angle. In this parametrization we get 

1-0(A 2 ) 0(A) 0(A 3 ) 

yCKM = I 0(A) l _ 0(A 2) 0(A 2) j (436) 

0(A 3 ) 0(A 2 ) 1-0(A 4 ) 
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Further, if we set 

P = 2.5 x w, 5 = 1.5 x w, (4.37) 

then we have quark masses 

m u = 0(X 7 v u ), m c = 0(X 3 - 5 v u ), m t = 0(v u ), (4.38) 
m d = 0(X 7 v d ), m s = 0(\ 6 v d ), m b = 0(\ 3 v d ). (4.39) 

These results are in line with the observed values. Since a, j3, 7, and 5 are set to be 
even positive integers, w should be a multiple of 4 in this case. Taking x w = A ~ 0.22 
and x 2sk ~ 2 = m 3 / 2 /M s = lCT^ 15 "" 16 ) into account, we obtain the constraint 

sk = (10 ~ 14) x w. (4.40) 

As a typical example, we will often refer the set of parameters 

s — 10, k — 5, w — 4, e = 0, 

a = 4, /3 = 10, 7 = 8, 5 = 6, (4.41) 

in which we have K = 51. In this case we have x ~ 0.7 and x k ~ 0.15. Consequently, 
when M s = 10 18 GeV, the symmetry breaking scales (S'o) and (Nq) turn out to be 
~ 7 x 10 17 GeV and ~ 1.5 x 10 17 GeV, respectively. Since the symmetry breaking 
scales are very large, we have the standard model gauge group over the wide energy 
range. 

5 Spectra of leptons 

Let us now study the mass matrices for lepton sector, in which L-H d mixing 
occurs at energies below the scale (Nq). Colorless S77(2)£-doublet fields L and H d 
are not distinguished with each other under G st . Then, hereafter we denote -R-parity 
odd H di as L\ (i = 1,2,3). As mentioned in section 2, H ui and L\ (i = 1,2,3) in 
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chiral multiplets do not develop their VEVs. It follows that there exist no mixings 
of SU(2)l x U(1)y gauge superfields with H ui and L\ [i = 1, 2, 3). Since both L and 
V are S , f/(2)^-doublets, the CKM matrix for lepton sector becomes a unit matrix 
irrespective of the magnitude of L-V mixing. For charged leptons the superpotential 
is 

W E ~ (SoS) h »S L' t H UJ + (S S) m » U(N c H uj + H d0 E*), (5.1) 
where the exponents are integers in the range < hij < K and satisfy 

2hij + hi + bj + a + 2 = mod K. (5.2) 

Thus we have 



2hij = 2h 33 + 



( 2(3 + 25 (3 + 25 (3 + 5 

(3 + 25 25 5 | mod isT (5.3) 
V (3 + 5 5 



with 2h 33 = —2b 3 — a — 2. As before, we introduce a 3 x 3 matrix H with elements 

H ij = 0{x 2h ^). (5.4) 

The mass matrix for charged leptons has the form 

H+ £ c + 

R V- ( xH \ (5-5) 
' L~ \ x k M p d M J 

in Ms units. This Mi is also a 6 x 6 matrix and can be diagonalized by a bi-unitary 
transformation as 

Vf'MM. (5.6) 



From Eqs. (|5.5| ) and (|5.6j ) the matrix 



ur'Mim a = ur 1 [ Al + B f g j n (5.7) 
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is diagonal, where 

Ai = x 2 H^H, Bi = r 2h .\I M. (5.8) 

The analysis is parallel to that of down-type quark masses in the previous section. 
We have the eigen equation 

det (A + B l - -^j = (5.9) 

for heavy states. For three light states their masses squared are given by the eigen 
equation 

det [x'^iAj 1 + Bf 1 )' 1 - = 0. (5.10) 

The light states correspond to observed charged leptons. Due to L-L' mixings mass 
pattern of charged leptons could be changed from that of up-type quarks. Introducing 
appropriate unitary matrices Wi and V«, we can diagonalize (A t + B{) and (Aj 1 + 
Br 1 )- 1 as 

Wi\A l + B^Wt = A; o) , VfW + BfY'Vt = A[ 2) , (5.11) 
where Ai -* and are diagonal 3x3 matrices. Masses of charged leptons are written 



as 



v 2 d (x-*AF>) u (< = 1,2,3). (5.12) 
Explicit forms of V t and Ui are 

^ " U fc MW,(A; o) )-V2 (x*M+)- 1 V I (A l (2) )V2 J' 

In the same way as the case of down-type quarks, we also choose a large mixing 
solution. Defining an even integer £ by 

£ = 2/i33-A; + l-(a + 7-/3-(5), (5.15) 
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we now impose the condition 

a + 7 + £-/3-5~0. (5.16) 

This condition means that 

2/i33~Jfe-l. (5.17) 

In this case eigenvalues of A[ + B\ become x 2<yk+l3+25 \ x 2 ( fc + 5 ) and x 2k . Thus masses 
of three heavy charged leptons are 

M s x k+I3+2S , .Us/''". M s x k , (5.18) 

whose eigenstates are mainly Li-H ul , L' x -H u i and (0(1)L' 2 + 0(1)L' 3 )-H u3 , respec- 
tively. On the other hand, light charged leptons have their masses of 

VdX *+P+y+s+^ v d x a+ ^\ v d x~<, (5.19) 

which correspond to observed e-, fi- and r- leptons, respectively. Main components 
of these eigenstates are L 2 -E\, L 3 -E 2 and (0(1) L' 2 + 0(1)L' 3 )-E 3 , respectively. The 
unitary matrix V/ which diagonalizes Aj 1 + Bf 1 , is expressed as 

/ i _ 0(x 213 ) 0(x?) 0(x f3+s ) \ 
Vi = 0(x 13 ) 1 - 0(x 25 ) 0(x s ) . (5.20) 
V 0(xP+ s ) 0(x 5 ) l-0(x 25 ) ) 

Taking 

i = w (5.21) 
together with the above parametrization ( 4.35Q and ( |4.37 ), we obtain 



m e = 0(X 8 v d ) J m M = 0(A 4 - 5 ^), m T = 0(\ 2 v d ) (5.22) 

and 



l-0(\ 5 ) 0(A 25 ) 0(A 4 ) \ 
V,= | 0(A 25 ) 1-0(A 3 ) 0(A L5 ) 

0(A 4 ) 0(A L5 ) 1-0(A 3 ) J 



(5.23) 
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We now proceed to study mass pattern of neutral sector. In the present framework 
there are fifteen neutral fields, i.e., H^, L[°, L°, Nf and Si (i = 1, 2, 3). For neutral 
fields the superpotential is of the form 

W N ~ (SoS^SoL'^ + (S S) m -L l (N^H UJ +H u0 N^ 

+ (SosyHs^iSjS) + (s sy*j(s t s)(NCN c ) 

+ (S S) n ^(N^N c )(N^N c ), (5.24) 
where the exponents s„, Uj and n»j are determined by 

2sij + a,i + dj + 2a + 2 = 0, 

2tij + (n + bj +a+b+2 = 0, mod if (5.25) 
2n„ + k + bj + 2b + 2 = 

in the range < Sij,tij,riij < K. From Eqs.( f4.2| ), ( f4.14| ) and ( |5.2| ) these equations 
are put into the form 

2sij = 2zij + a - 2a, 

2tij = 2rriij + b — a — b, mod K (5.26) 
2riij = 2hij + ao — 2b. 

In the above we have imposed the ansatzs 

2m 33 = 0, (5.27) 

2z 33 = k-l-a--y + P + 8, (5.28) 

2/i 33 = k-l + a + j-(3-5 + £~k-l (5.29) 

together with the conditions 

a + a = 2, b + b = 2k, a + 2b = 2e, e = 0, 1. (5.30) 
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Therefore, the exponents sy, Uj and rewritten as 



2sjj — 2/2<2jf 2 "I - 2c 2/c 

2t y = 2mij - 2 + 2e - 2fc, mod If 

2n»j = 2/ijj + 2e - 4A;. 

Introducing 3x3 matrices S, T and N with elements 



(5.31) 



(5.32) 



we have a 15 x 15 mass matrix 



HI L 



10 



a 



M 



N 



N c S 

\ 



Pu M 

p u M T x 2k N x k+1 T T 

k + lrp Jl 



El ( xH x k M T 
L'° xH 
L° x k M 

N c 
S 



(5.33) 



\ x k+l T x 2 S J 
in M s units for neutral sector, where p u = v u /M s . Since SU(2) L symmetry is pre- 
served above the electroweak scale, the eigen equation for six heavy states is the same 
as Eq. ( |5.9| ). For another nine states we have an approximate eigen equation 

V 



det M LNS - 



0. 



(5.34) 



where Mlns is defined by 



M, 



LNS 



( 



\ 



e^AfV^Vr 1 \ 

s„V,(A[ 2) ) 1/2 x 2k N x k+1 T T 
x k+1 T x 2 S ) 



(5.35) 



with e u = x k p u . Light neutrino masses are given by 



m. 



M s x 2k \VdJ v Jn 



(5.36) 



with 



A N = (N — T T S^Ty 1 . 



(5.37) 
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This type of mass matrix has been discussed in Ref . ]TJ| . Let us suppose that M 



iLNS 



is a 3 x 3 matrix. When N — TS l T ~ A" in order of magnitude, the usual seesaw 



mechanism |l7j is at work. On the other hand, when iV < TS~ 1 T, another type of 
seesaw mechanism takes place. In the present framework Mlns is a 9 x 9 matrix. 
The exponents 2n,ij in are equal to those in (T T S~ 1 T) i j in modulus K. The 
relative magnitude of iVy to (T T 5' _1 T)jj depends on the value of k. In the case 
k > §(a: + 7) — 1 + e which corresponds to small s (s 8), light neutrino masses are 
not so small. For instance, we have m VT = O(lkeV). So we do not adopt this case. 
In the case < |(a + 7) — 1 + e which corresponds to s <: 8, light neutrino masses 
are extremely small. Specifically, the latter case yields 

fVf'AjvV,) = J x 3fc+l-2(/3+5 + e) ; x 3/c + l-2( 5+ e) ; a .3fc+l-2el _ (g gg) 

V / ii 

Combining this with the result ( |5.19|) for ran, we have the light neutrino masses 



mui = ^x k+1 ~ 2e x {x 2 ^ a+ ^\ x 2 ^\ x 2 - 1 }. (5.39) 
M s 



The previous parametrization ( |4.35|) and (|4.37| ) leads us to 



It follows that 



m vi = x k+1 ~ 2e x {A 8 , A 6 , A 4 }- (5.40) 

M s 



m - ^ TZE „ mvT < O(10" 7 eV). (5.41) 



A 4 A 2 

The calculated neutrino masses seem to be too small. According to the analyses 



of solar neutrino fl8 |, atmospheric neutrino [|l^] and cosmological constraints |20|, it 
is preferable that three typical mass scales of neutrinos are ~ 10~ 3 eV, ~ 10 _1 eV 
and ~ lOeV. The ratios m ve jm VVL and m UfJi /m UT obtained here are consistent with 
those among the above three typical mass scales. As pointed out at the beginning 
of this section, the CKM-matrix for lepton sector is a unit matrix irrespective of 
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the magnitude of L-V mixing. The situation is unchanged even through seesaw 
mechanism. This is because the fields N c and S are S'C/(2)x,-singlet. In addition, the 
components of N c and S in light neutrinos are very small. Thus we have no flavor- 
changing charged currents at tree level. Recently, by introducing discrete symmetries 
which yield appropriate texture zeros in Yukawa couplings, the hierarchical pattern of 
neutrino masses has been examined in Ref . . Finally, we touch upon the remaining 
eigenvalues of Eq. (|5.34| ). Three pairs of heavy states which are G st — neutral have their 



masses of 

M s x~ k - 1+2e x {x p+6 - a -t, x a+ \ x a+7+5+f }. (5.42) 

6 Proton decay 

After studying the particle spectra of down-type colored fields in $ an d 
we explore the proton stability in this section. Under the discrete .R-symmetry the 
superpotential of down-type colored fields is of the form 

W g ~ (S SyS g g c + (S Sy k -^(SgT + ( W'^o^gr) 

+ (S Sy k ~ 1 (g g + g c o r + (So^'WW) + (S S) {s - 2)k+e - 

+ (S S) {s - 1)k D^D c + (S SY s - 2)k+q+e+1 SoN c g^D c , (6.1) 

where q = k + |£ — e — 2. When So, S, and N° develop non-zero VEVs, the 
superpotential induces the mass matrix 

go ( 0(x 2k - 3 - 2e+ t) 0{p) 0(px k+2 ~ 2e ) \ (6 2) 

0( px i-fc+C) 0(px~ 3k+2e ) 0{px- 2k+2 ) 
g c \ O(p) 0(px- 2k - 1+2e -t) 0(px- k+1 -Z) 



M g = D 



in Ms units, where p = m^^/Ms = 10 ( 15 ~ 16 ). This mass matrix is diagonalized by 
a bi-unitary transformation as 

V g l M 9 U 9 (6.3) 
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with 



I x 0{px-' ik+A+2e ) 0(px- 2k+3+2e ^) 

0(px- 3k+4+2e ) 1 ©(ar*- 1 -^) | . (6.4) 

V 0(/)x- 2fc+3+2e "«) o{x k - 1 -^) 1 

/ 1 0(par 2fc+3+2e ) 0(px- k+5 -t) 

0(px- 2k+3+2e ) 1 0(x fe+2 - 2e ) | . (6.5) 

V 0(px- k+5 S) 0(x k+2 - 2e ) 1 



The eigenvalues are given by 



Mg A 


= 0(M s x 2k 




Mg B 


= 0(m 3/2 x~ 




M g c 


= 0(m 3/2 x~ 





(6.6) 



From explicit forms of V g and U g we find that three eigenstates gA~g% 9b-9 c b an d 
9c~9c are approximately go-g^, D°-g and ~g c -D§ states, respectively. In the typical 
example ( |4.41| ) M 9 a is nearly GUT-scale (~ 10 16 GeV). By contrast, M 9 b and M g c are 
as small as O(10 5 ~ 6 GeV). Then, at first sight, it seems that dimension-five operators 
mediated by these rather light colored fields lead to fast proton decay. However, this 
is not the case. The dimension-five operators mediated by light colored fields are 
strongly suppressed because of extremely small effective couplings. This is due to the 
fact that (1,2), (2,1), (1,3) and (3,1) entries of V g are O(10 -13 ). In what follows we 
explain this situation more explicitly. 

Since i?-parity of quark and lepton superfields are odd, those of S*?7(2) L -singlet 
colored superfields mediating proton decay should be even. The relevant superfields 
are g , g£, and Dq which reside in <3> (27). Effective trilinear couplings with g , g§ 
and -Dq are given by 

Wg ff ~ (Q T ZQ + N cT HD c + E cT HU c )g + (Q T ZL + U cT HD c )g c 

+ (Q T ML' + D' T MN C + D' cT MU c )D c , (6.7) 
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where 3x3 matrices Z, H and M have already been determined in sections 4 and 5. 
Superfields go, g$ and Dq are expressed in terms of mass eigenstates as 

90 = Vgll 9 A + Vgl2 9B + Vgl3 9C, 

9o = Ugug C A+Ugi2 9 C B+Kg l3 g c c , (6.8) 
D c = Ug 31 g c A +U g32 g c B +Ug 33 g c c . 

Taking V g u, U g \\ — 1 an d U g3 \ ~ px~ k+5 ~^ into account, we can obtain dominant 
dimension-five operators from gA-g% exchange 

_L x -2k + 3 + 2e~^ Q T ZQ + N cT HD c + £] cT ffJJ c ^Q T Z L + U cT HD°). (6.9) 

M s 

Similarly, dominant dimension-five operators from gc-g c c exchange are 

— x - k+2+2e (Q T ZQ + N cT HD c + E cT HU c ) (Q T ML' + D' T MN C + D' cT MU c ). (6.10) 
M s 

The prefactor is induced from V ff i3(M g cO -1 W ff 33. In dimension-five operators for gB-g% 
exchange the prefactor is V g i 2 (Mg B )~ 1 Ug 32 , which is smaller than the one for gc-g°c 
exchange by the factor x 2fc + 4 ( 1 - e ). Therefore, the study of dimension-five operators 
coming from gA-g°A an d Qc-Qc exchanges suffices to explore the proton stability. 

We now rewrite the above operators in terms of quark and lepton mass eigen- 
states, which are represented by using the symbol "tilde". In order to implement this 
rewriting, we can use the transfer 

Q -> V U Q' =( V q y A I yCKMf) ) > 



' u 

U c -> U U U C , D c ^U d22 D^ 

D' - V d D, D' c ^U dl2 Dc, (6.11) 
L — > Vi 22 L, U — > Vn 2 L, 



E c -> V/£ c , 



27 



where 



UdYl — 


-(xZ)- 1 V,(Af ) ) 1/2 , 


(6.12) 


U<m = 




(6.13) 


Vll2 = 


-(^t)-V,(Af))V2 ) 


(6.14) 


Vl22 = 


(x fc Mt)- 1 V i (Ap ) ) 1/2 , 


(6.15) 



which are three- by-three blocks of matrices given in Eqs.( |4.24"D and ( |5.13| ). Light 



component of A^ c is extremely small and then its contribution to nucleon decay is 
negligible. Therefore, Eq. (|6.9| ) is translated into 



±- x -^^[QV{V T u ZV u )Q' x Q' T (V^ZV l22 )L 

M s 

+ E cT (VfHU u )U c x U cT (U^HU d22 )D c }. (6.16) 
Similarly, Eq. ( |6.10| ) is put into 

±- x -^[QV{VlZV u )Q> x Q' T iy T u MV ll2 )L 

M s 

+ E cT (VfHU u )U c x U cT (UlM T U dl2 )b c }. (6.17) 
The dimension-five operators result in nucleon decay via gaugino- or Higgsino-dressing 



processes [p2| ||23|| . Among various dressing processes the exchange of charged wino 
or Higgsino give predominant contributions to nucleon decay. Since S77(2)£-singlet 
states U c , D c , E c do not couple to S'f/(2)i-gauginos, dominant dimension-five opera- 



tors with charged wino-dressing processes turn out to be the first terms in Eqs. ( |6.16| ) 



and ( |6.17| ). Thus we have dominant operators incorporating charged wino-dressing 



processes 



1 x - 2fc+ 3 + 2e- 5 g Uv T ZK)ii g ; x Q> (yj^^ 
Ms 

+^r x ~ k+2+2e Qi^uZv u ) ll Q' 1 x Q' 2 (ylMv ll2 ) 2] L r (e.is) 

M s 
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Simple calculations yield 



■I J 



0(x k - 1+ t)(VTMV ll2 h 

0(x a+ ^) 0(x a+ ^) 0(x~<) 

x k-l+l3+8 x I ( X 7+?) O(X^) 0(x~<- a ) 

0(x*) 0(1) 0(x~ a ) 



(6.19) 



(6.20) 



From these relations Eg. ( |6.18| ) becomes 
1 



Mc 



(Q' i g' i g' 2 L,)x 2 ( Q+ ' 3+ ^ +1+2e -« x (x- a+ Z, x ~ x 



ct ^—2a\ 



(6.21) 



This implies that a dominant mode of proton decay is p — > K° + In this decay 
mode the magnitude of the dimension-five operator is given by 



J_ a+2(p+i+8)+l+2e-£ ^ _L,\ 12 > 10" (25 ' 5 ~ 26 ' 5) GeV _1 

M s M s 



(6.22) 



The second terms in Eqs. (|S.16|) and (|6.17| ) contribute to nucleon decay via charged 
Higgsino-dressing processes. The relevant terms are 
1 



-x 



+ 



Mc 



-X 



-k+2+2e pc/i ;T 



(6.23) 



where (J, m) = (1, 2), (2, 1). Using the relations 



(Idl HUd22)ij 



0{H„ 



Oix^XU'MU^h 

( 0(x a+l3+s ) 0(x a+l3+s ) 0{xP +s ) 

0(x a+s ) 0(x a+s ) 0(x 5 ) 

V 0(x a ) 0(x a ) 0(1) 



x k - l+ ~< X 



(6.24) 



(6.25) 



we obtain the dimension-five operators for 5 , f/(2)^-singlet fields 
1 



-{E c MU0l) x p+2{a+ ^ +5)+1+2e x (1, x~P, 



(6.26) 
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In the Higgsino-dressing processes the operators are multiplied by their Yukawa cou- 
plings. As a consequence we have a dominant contribution in the case j = 3. The 
above operator multiplied by the Yukawa couplings for E% and C/f becomes 



in magnitude and results in the decay p — > K + + V T . 

In conclusion of this section, the main mode of proton decay is p — > K° + 
in which the magnitude of the dimension-five operator is about \ 12 /M S ~ 
1 Q-(25.5~26.5)Q ev -i > rp^^ j m pij es triat trie p ro ton lifetime is about 10 33 ~ 35 yr. This 

result is consistent with the present experimental data. 

7 Gauge coupling unification 

As is well-known, there is a discrepancy between the string scale Ms and the 
MSSM unification scale ~ 2 x 10 16 GeV. Main concern here is whether or not we 
can reconcile this discrepancy in the present model. For this purpose we study the 
renormalization group evolution of the gauge couplings in the model up to two-loop 
order. 

In the preceding sections particle spectra have been already studied. Unlike the 
MSSM, in the present model there are many extra intermediate-scale fields, which 
are tabulated in Table II. In particular, the contributions of (H u q, Hdo, H u , Hd) and 
(<7, ~g c , D c , Dq) are significant, because their masses are lying in rather low energy 



The evolution equations for ctj = gf/'in are generally given up to two-loop order 




.2(a+2 7 +5)+l+2e 




(6.27) 



region (10 2 ~ 7 GeV). 



Table II 
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by 



dfi 2 



1 

All 



47T~ 



47T 



(7.1) 



where fj, is the running mass scale [ 24jl . The coefficients bi, by and Oj are determined 
by the particle content of the model. The third term in the r.h.s. represents the con- 
tribution of Yukawa couplings. In the present calculation we take accout only of the 
largest Yukawa couplings / = M 33 , namely, Yukawa couplings of the third generation 
$3(15, 1)^3(6*, 2)$o(6*, 2) and for simplicity we neglect the renormalization group 
evolution of the Yukawa couplings. In our analysis it is assumed that string threshold 
corrections are negligibly small. 

In the region between M5 and (So) = Ms x, where the gauge symmetry is SU (6) x 
SU (2)r, we have 



9 



9 15 
175 81 



28 
60 



(7.2) 



where y = f 2 /Air is taken to be a constant. In the region between M$x and (Nq) 
M$x k , where the gauge symmetry is SU(4)ps x SU(2)l x SU(2)r, we get 









( 118 


9 


15 






bij = 


45 


53 


15 








V 75 


15 


81 



0£ 

y 



36 

32 + An(S 3 ) 
56 + An(S 3 ) 



(7.3) 



with 



n(S 3 



1 M s x> n>M s x - k+ ^ +2e - 1 
M s ar fe +7+2e-i > ^ > M; 



sx k . 



(7.4) 



In the wide energy region ranging from Msx k to m 3 /2 = Ms x 2sk ~ 2 the gauge group 
coincides with the standard model gauge group. From Table II we can calculate the 
coefficients, which are of the forms 







J I 




- 






I 6 / 


V 2/5 




(7.5) 
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J K1 



/ 14 9 11/5 

24 18 6/5 

V 88/5 18/5 38/5 

/ 



+ 



+ 







7 3/5 
V 9/5 9/25 



/ 34/3 4/15 


V 32/15 8/75 



(7.6) 



y 





f 4 ^ 




8 ^ 




6 


n{Dl) + 


3 n(D c Q ) 


v 44/5 ; 


V 14/5 J 




{ 31/5 ; 




n(iV 3 c ). 



(7.7) 



In these expressions nn and n 9 stand for the numbers of doublet Higgses and extra 
down-type colored fields, respectively and are given by 



n H = < 



= < 



4 
3 
2 
1 

5 
4 
3 
2 
1 




n(D^), ti(Dq) and n(N£) are 
n(D 3 c ) = | 
n(D c ) = | 
n(iVJ) = { 



k+/3+8-a 



M S X k > fl> M S x k+s 

M s x k+s > /i > M 5 x fc+/3+M 

M5X fc+/3+25 > ^ > MsX 2sk-Ak+2e-l 

M s x^-^+ie- 1 > ^ > M 5 x 2sfc - 2 , 

M S X k > fJL> M S a;2fc-3-2e+5 
MsX 2k-3-2eH > jj>M S X 

M s x k+li+s - a > fx>M s x k+a+ ^ +s 
M s x k+a+ ^ +s > fi>M s x 2sk ~ 3k+2e - 2 
M s x 2sk - 3k+2e - 2 > fi> M s x 2sk - k - 1 -t 

MsX 2sk-k-l-t > ^ > M 8 X 2sk ~ 2 . 



M s x k > n> M s x k+l3+5 - a 
M s x k+ P +5 ~ a > n > M s x 2sk - 2 , 

M s x k > n> M s x 2sk - k ~ 1+i 
M s X 2sk ~ k - 1H > fj, > M s x 2sk ~ 2 , 

M s x k > /2 > M s x a+r, - k+2e - 1 

MsX a +1 -k+2e-l > ^ > MsX 2sk-2_ 



(7.8) 



(7.9) 



(7.10) 
(7.11) 
(7.12) 



It should be noted that in the present model we obtain nn — n g = over rather 
wide energy range. By contrast, in the MSSM we have n H — n g — 1. In the region 
between m 3 / 2 and M z where supersymmetry is broken, all superparticles except for 
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light Higgses do not contribute the evolution equations. This leads us to 





( - 7 \ 


/ -26 


9/2 




-7/3 , 


K = 12 


97/6 




I 23/5 / 


V 44/5 


9/2 



y 




We are now in a position to solve the evolution equation numerically. The behavior 
of the renormalization group flow is shown in Fig.l, in which we choose a typical 
example ( [4.41 ). In the present calculation the parameters are taken as 



M s = 0.5 x 10 18 GeV, m 3/2 = 200GeV, aj t l mg = 14.0, / = 1.7, (7.14) 

where a str ing represents the unified gauge coupling at the string scale. Resulting 
values of a~ 1 {M z ) are 



a^(Mz) = 58.96, a^[M z ) = 26.03, a^(M z ) = 8.66. 



(7.15) 



Compared with the present experimental values |25 



a^(M z ) = 58.95±0.08, a^(M z ) = 29.66±0.07, a^(M z ) = 8.48±0.43, (7.16) 

the calculated a 1 " 1 (M^) and a^ l {M z ) are consistent with the data, while the cal- 
culated oo 2 x [M z ) is smaller than the observed one by ~ 3.5. Consequently, our 
analysis shows that the gauge coupling unification is not achieved at the string scale. 
This suggests that we are not successful in getting proper particle spectra of extra 
intermediate-scale matter. 



Fig. 1 
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8 Summary 



In the context of level-one string model we have explored a possibility that 
characteristic pattern of quark-lepton masses and the CKM matrix have their origin 
in the discrete /^-symmetry and mixing mechanism. In this paper we have chosen 
Zk x Z 2 symmetry with K = sk + 1 as an example of the discrete i?-symmetry. The 
-^-symmetry is assumed so as to be in accord with the i?-parity in the MSSM and 
is unbroken down to the electroweak scale. The vector-like multiplets $ , $ and the 
chiral multiplets $j (i = 1,2,3) are assigned to even and odd i?-parity, respectively. 
Under this assignment no mixing occurs between the vector-like multiplets and the 
chiral multiplets. The Z K symmetry is used as a horizontal symmetry. The Z K 
symmetry controls a large hierarchy of the energy scales of the symmetry breaking and 
particle spectra. Triplet-doublet splitting problem and the /z-problem are solved as a 
result of the discrete symmetry. The assignment of -/^-charges to chiral multiplets is 
of great importance in explaining the observed hierarchical pattern of quark-lepton 
masses. 

The mass hierarchy of up-type quarks is a direct result of the horizontal discrete 
symmetry. On the other hand, for down-type quarks there appears a mixing between 
D c and D' c {= g c ). Due to the maximal mixing mass pattern of down-type quarks is 
different from that of up-type quarks. The mass hierarchy obtained here is 

m u = 0(v u x a+ ^ +s ), m c = 0{v u x^ 6 ), m t = 0(v u ), (8.1) 
m d = 0(v d x a+ ^ +5 ), m s = 0(v d x^ +5 ), m b = 0(v d x~ a+ ? +5 ). (8.2) 

These results are consistent with observations under the parametrization a — w, (5 — 
2.5u>, 7 = 2w, 5 = 1.5w and x w = A ~ 0.22. Further we obtain a phenomenologically 
viable CKM matrix. In lepton sector L-L'(— H d ) mixing occurs. Hence, mass pattern 
of charged leptons is also changed from that of up-type quarks under a large mixing. 
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The obtained mass hierachy for leptons is 



m e = 



0(v d x a+ ^ +s+ t), 



0(v d x a+ i +s ), 



m T = 0(v d x J ). 



(8.3) 



The CKM matrix in lepton sector amounts to a unit matrix irrespectively of the 
magnitude of L-L'(= H d ) mixing. This is because both L and L'{— H d ) are SU(2) L - 
doublets. Therefore, lepton flavor violating processes are extremely suppressed. See- 
saw mechanism is at work for neutrinos. For large s (s > 8) light neutrino masses 
are 



These masses seem to be too small compared with those expected from solar neu- 
trino and atmospheric neutrino data. In the present framework the proton lifetime 
is 10 33 ~ 35 yr, which is long enough to be consistent with experimental data. The sup- 
pression of the dimension-five operators occurs because of the superheavy mass of the 
mediating particle for certain processes and because of the extremely small couplings 
for the other processes. On the other hand, we are not succesfull in achieving the 
unification of gauge couplings at the string scale. Nevertheless, it is suggestive that 
the obtained numerical value a~ t l ing ~ 14 corresponds nearly to the self-dual point 
g s tring = 1 with respect to 5-duality (strong/weak duality). 

Both in D c -D' c (= g c ) and L-L'(= H d ) mixings the mass differences between heavy 
states and light states are extremely large in order of magnitudes. This implies 
that these mixings do not practically bring about flavor- changing neutral current 
processes. In addition, flavor-changing neutral current processes via superparticle 
exchanges at loop level are also suppressed enough to be consistent with experimental 



m 



m 



m 
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data, provided that the soft SUSY breaking parameters are universal at the string 
scale. More explicitly, the most stringent experimental bound on the mass difference 
of squarks d and s is derived from the K°-K° mixing. As pointed out in section 
4, S77(2)i,-singlet components of down-type quarks d, s and b are nearly 0(1)D\ + 
0(1) D%, D'i and _D 3 C , respectively. Although D c and D' c are indistinguishable from 
each other under the standard model gauge group, D c and D' c reside in (6*, 2) and 
(15, 1) of SU(6) x SU(2) R , respectively. Further, D c and D' c reside in (4*, 1, 2) and 
(6, 1, 1) of SU(4:)ps x SU(2)l x SU(2) r , respectively. Therefore, gauge interactions 
cause soft SUSY breaking masses of d c R and s R to evolve differently through radiative 
corrections in the energy region ranging from M$ to M$x k . However, in the present 
model this energy range is rather nallow. In fact, x k is about 10 _a8 in a typical 
example. Consequently, it can be shown that the difference Sm 2 = m 2 (d c R ) — m 2 (s c R ) 
remains small at low energies. Let us estimate numerically the difference 6m 2 in a 
typical example. When we assume 5m 2 (Ms) = 0, the difference at the scale M$x 
(= (So)) becomes 

5m 2 (M s x) ~ -0.016 x M\ (8.5) 

through the RG evolution, where Ma is an averaged gaugino mass. Subsequently, 
the RG evolution from M s x to M s x k (= (N£)) leads to 

5m 2 (M s x k ) - 5m 2 (M s x) ~ 0.008 x M\. (8.6) 

Combining these two results, we obtain 

5m 2 (M s x k ) ~ -0.008 x M 2 A . (8.7) 

Since Yukawa couplings of down- type quarks are tiny in case of tan /3 ~ 1, the contri- 
butions of Yukawa interactions to Sm 2 (7713/2) are small compared with Sm 2 (Ms x k ). It 
follows that Sm 2 (m 3 /2) — 5m 2 (M s x k ), which is consistent with a bound on Sm 2 (7773/2) 
given in Ref . |26[ . 
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Although we did not deal with CP-violation, there are two possibilities of intro- 
ducing the CP-phase in the present framework. One possibility is that the CP-phase 
comes from complex VEVs of moduli fields. In this case the coefficients of the terms 
in the string-scale superpotential are complex in general. Another possibility is the 
case that the coefficients in the superpotential are all real but VEVs (So), (S), (JVq) 
and (N c ), are complex. When we take the relative phase of (Nq)(N c ) to (So)(S) into 
account, there appears CP- violating phase in the model. 
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Appendix A 

In this appendix we show that the minimization of the scalar potential yields tree- 
level breaking of the gauge symmetry under an appropriate condition on soft SUSY 
breaking parameters. In minimal supergravity model the soft SUSY breaking terms 



are given by 27 



ft 



m 3/2 9 2 B + m 3/2 OB*- m\ j2 2 9 C exp(2gV)$ 



J d 2 6m 3/2 6 2 AW + h. c. 



(A.l) 



Here 7713/2 is supposed to be O(lTeV). The universal soft SUSY breaking parameters 
A, B and C are generally zero or order unity. Although A and B are generally 
complex numbers, C is a real one. This type of £ so ft leads to the scalar potential 

2 

+ m 3/2 {AW + A* W*) 



V 



E 



dW 



+ m 3/2 J2(B(f ) ~ + B* 



+ m 2 3/2 (C + \B\ 2 ) J2 \U 2 + p-tenn), 



(A.2) 



where fa is a scalar component of the chiral superfield $j. In the above expression 
it is assumed that the terms of higher powers of 1/M p i are negligibly small. We will 
shortly show that this assumption is justified. 

For illustration, we take one set of vector-like multiplet $ and <3>, whose scalar 
components are denoted as and 0, respectively. Let us consider the case that the 
nonrenormalizable interaction 



W = AM|- 2n ($$) n 



(A.3) 



is compatible with the discrete symmetry, where A is a positive 0(1) constant and 
n is a large positive integer. In a typical example (|4.41|) we put n = sk = 50. For 
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simplicity we denote dimensinless quantities V/Mg and 4>/M s by the same letters as 
the original V and 0. Thus 



V 



^2 \ 2 

n A 



in-r 



+ 



+pA {(A + 2n£?)(00) n + (A* + 2nB*){<p*f) n ) 
+p 2 (C + |5| 2 )(|0| 2 + |0| 2 ) + (D-term) 



(A.4) 



with p = m 3 / 2 /M s . Minimization of V leads to the D-flat direction 



= x. 



Writing the phase factor of VEVs explicitly as 



(A.5) 



(0)(0)=xV e , 



(A.6) 



we have the scalar potential 



V = 2n 2 Yx An - 2 + 2p\A + 2nB\ \x 2n cos5 + 2p 2 (C + \B\ 2 )x 2 , (A. 7) 

where 5 = n6 + arg(A + 2nB). From the stationary condition dV/dS = and x > 
the phase 6 is determined as cos S = —1. Therefore, the dependence of V on x is 
given by 



V = 2 



n\x 



2n-l 



B + — 

2n 



2p 2 



B+ — 

2n 



C- 151 



x 



(A.8) 



Consequently, if the inequality 



B + — 

2n 



>C+\B\ 



(A.9) 



holds, V is minimized at a nonzero value of x, namely at x ~ p 1 /( 2n ^ 2 ). If C < 0, 
the above inequality is satisfied, for example, in the case |arg(A_B*)| < ir/2 even 
for large n. It is worth emphasizing that the soft SUSY breaking mass parameter 
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(C + \B\ 2 ) is not necessarily negative. If only the above inequality is satisfied, the 
gauge symmtry is spontaneously broken at tree level. It is not necessary for us to 
rely on the radiative symmetry breaking mechanism. In this paper the exponent n 
is taken to be rather large. The larger n implies the larger VEV \{(j>)\ = M s x. The 
large value of | (0) | is consistent with the tree-level symmtry breaking. 

In supergravity theory with canonical Kahler potential the supersymmetric term 
of the scalar potential is expressed as 

2 



V = e 



K/M 2 



E 



d(j) t M 2 



M 2 



\W\ 



+ (D— term) 



(A.10) 



with M = M p ;/V87r ^ M s . In the present model we get 



i 8 —) 

v 84, 1 



~ n\x 2n ^M 2 , 



\x 2 ^)Mi 



(A.ll) 



Since n is large and x < 1, M s /M ^ 1, V is dominated by dW/dcf). The overall 
factor (exp(K/M 2 )) is order unity. Therefore, the above analysis is relevant to the 
issue of the symmetry breaking. 
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Appendix B 



In this appendix we address to the issue of -R-parity conservation within the present 

framework. It is shown that if eigenvalues of the mass matrix 

' x 2k N x k+l T T \ 
x k+1 T T x 2 S 



M 



NS 



(B.l) 



in Ms units are sufficiently large compared with m 3 / 2 , the scalar potential is mini- 
mized along the direction where -R-parity is conserved. The mass matrix M NS is a 
submatrix of Mlns given in section 5 and yields masses of -R-parity odd and G st - 
neutral superfields. It has already been found in section 5 that the above condition 
is satisfied for the solutions discussed in the text. 
The superpotential can be separated as 



w = w 1 + w 2 , 



(B.2) 



where W 1 is a function only of -R-parity even fields So, S, N$ and ISF, while each 
term of W 2 contains -R-parity odd fields $j = Sj, N? (i — 1, • • • , 6 ; j — 1, 2, 3). In 
the same manner as the notations in appendix A, we now use dimensionless quatities 
in M s units. Due to the Z K x Z 2 symmetry the explicit form of W\ is given by 



W^^CriSoSf-^iN^Ny, 



(B.3) 



r=0 



where c r are 0(1) constants in Ms units. This superpotential satisfies a relation 



W l 



2sk 



So^ + S W 



dS 



dS 



+ 



2s 



8Wa —cdWi 
m 4- n 



W 2 is a even function of <£>;. Consequently, the scalar potential is of the form 



V = V l + v 2 



(B.4) 



(B.5) 



with 



V 1 



+ 



(So ^ S) 



41 



+ 



dW A 



dNS 



+ 



V2 = E 



i=l 



~ P 
-P 

dW 2 



B + 
B + 

YpB*^> 



A 



2sk 
A 2 



2s 

2 



-C-\B\ 2 )(\S \ 2 + \S\ 2 
C-\B\ 2 ^\N^\ 2 + \N C \ 2 ^, 



+ p(AW 2 + A*W 2 *)+p 2 Cj2\& 



(B.6) 
(B.7) 



i=l 



Here 0j's (i — 1, • • • , 6) represent scalar components of Sj and N? (j = 1, 2, 3) and 
W 2 is defined by 



W 2 = W 2 



2sk 



dw 2 ls dw 2 



1 

2s 



U a»rc 



dN§ 



dN 



(B.8) 



Scalar components of So, S, Nq and N are denoted by the same letters as the 
superfields themselves. As discussed in appendix A, under the assumption 



B + 



A 



2sk 



B + — 

2s 



>C+\B\ 2 , 



(B.9) 



So, S, Nq and N° develop nonzero VEVs and then the gauge symmetry is sponta- 
neously broken at tree level. The stationary condition is satisfied at nonzero values 
of S'o, S, Nq and N° and vanishing At this stationary point we get a negative 
value of the scalar potential 



V = V 1 = -O(p 2 (S > 2 ). 



(B.10) 



The question here is whether this point is the absolute minimum or not. 

Let us suppose that some of 4>i develop nonzero VEVs at the absolute minimum 
point. For such (pi, if 

»p|(0,)l, (B.H) 



then V 2 is dominated as 



( &W 2) 



(B.12) 
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and lifts up the scalar potential V. It follows that this point can not be the absolute 
minimun. Therefore, the relation 



, dW 2 



Pi <</>;> I 



(B.13) 



should be satisfied for all i. On the other hand, the mass matrix of </>j(<J>j) is given by 



d 2 Wo 



■) = (m ns ) 



(B.14) 



This matrix yields masses of -R-parity odd and G^t-neutral superfields, which are 
assumed to be sufficiently larger than p = m 3/ / 2 /Ms. Namely, when we introduce a 
unitary matrix Uns which diagonalizes M/vs, this assumption is expressed as 

d 2 Wo 



-) (u NS ) ki »p 



(B.15) 



for all i. Although we have six unknown parameters (4>i), there are twelve constraints 
( |B.13| ) and (|B.15|) on (0,) in which the orders of magnitude are quite different. Since 
we have too much constraints on (<fii), in generic case there are no consistent solutions 
except for = for all i. Consequently the absolute minimum of V is achieved at 
(<pi) = 0. This means that -R-parity is conserved. 
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Table I 





$ $ 


$1 $2 $3 


(15,1) 
(6*, 2) 


(a, +) (ao, +) 
(b,+) (b ,+) 


(ai,-) (a 2 ,-) (a 3 ,-) 
(&!,-) (6 2 ,-) (6s,-) 



Table II 



i?-parity 


Matter fields 


X (mass scale: m = 0(M s x x ) ) 






1 




f/ c , e c , r, r, 


k 




7^0 + 5) 


2sk - 2 


+ 


72(^0+ AO 


2sk - 2k 




-f^uO, -f^dO, H u , Hd, 


2sk-Ak + 2e- 1, 
2sk - 2e + 1 




9o, 9o, D%, g, g c , D° 


2k - 2e - 3, 
2sk - 3k + 2e - 2, 
2sk — k — l—£ 




H u i, L\ 


k + 13 + 25, k + 5, k 




D[, D'f 


k + a + 7 + 5, k + (3 + 5 — a, k 






-fc + 7 + 2e- 1, 
-fc + a + 7 + 2e- 1, 
-fc + /3 + 25 + 2e - 1 




Qi, U9, D9, Li, E\ 


> 2sk - 2 
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Table Captions 



Table I The numbers and bi (i = 0,1,2,3) in the parentheses represent 
the Zfc-charges of chiral superfields $(15, 1) and $(6*, 2), respectively, a and b stand 
for those of mirror chiral superfields $(15*, 1) and $(6,2), respectively. Respective 
Z 2 -charges (.R-parity) of the superfields are also listed. 

Table II Particle spectra in the present model. The number X stands for 
the exponent of x for the mass scale m = 0(M s x x ) of each superfield. Note that 
x 2sk-2 _ m3 / 2 /Ms and K — sk + 1. The parameters a, (3, 7 and 5 are given in 
section 4. In this table Df and D'f (Li and Z^) stand for light and heavy eigenstates, 
respectively, which are derived via mixings between D\ and D'f (Li and Z^). 



Figure Captions 

Fig.l The renormalization group flow of gauge couplings. The string scale, 
the soft SUSY breaking scale and the unified gauge coupling are taken as M s = 
0.5 x 10 18 GeV, m 3 / 2 = 200GeV and a~ t l ing = 14.0, respectively. The Yukawa coupling 
of the third generation is fixed to / = 1.7. 
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